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Abstract 

Finite-dimensional Lie algebras of vector fields determine geomet- 
rical collective models in quantum and classical physics. Every set of 
vector fields on Euclidean space that generates the Lie algebra sl(3,R) 
and contains the angular momentum algebra so(3) is determined. The 
subset of divergence-free sl(3,R) vector fields is proven to be indexed 
by a real number A. The A = solution is the linear representation 
that corresponds to the Riemann ellipsoidal model. The nonlinear 
group action on Euclidean space transforms a certain family of de- 
formed droplets among themselves. For positive A, the droplets have 
a neck that becomes more pronounced as A increases; for negative A, 
the droplets contain a spherical bubble of radius |A|3. The nonlinear 
vector field algebra is extended to the nonlinear general collective mo- 
tion algebra gcm(3) which includes the inertia tensor. The quantum 
algebraic models of nonlinear nuclear collective motion are given by ir- 
reducible unitary representations of the nonlinear gcm(3) Lie algebra. 
These representations model fissioning isotopes (A > 0) and bubble 
and two-fluid nuclei (A < 0). 
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1 Introduction 



The geometrical models of nuclear rotations and vibrations are based upon 
the general linear group GL(3,R), the group of invertible linear transforma- 
tions of three dimensional Euclidean space, and its subgroups, the special 
linear group SL(3,R) and the rotation group SO (3) [|l], H, If. These three 
groups are the motion groups of their corresponding models, the general 
collective motion model GCM(3) 0, ^, the special collective motion model 
SCM(3) [g 0], and the adiabatic rotational model R0T(3) §, ||. 

A collective model describes a dynamical system if the particle trajectories 
are compatible with the model's motion group. For example, the quantum 
adiabatic rotational model and the classical Euler rigid body model are useful 
when the velocity v of a particle located at the position vector x is given by 
V = uj X X, where u is the angular velocity of the rigidly rotating system. In 
this case we say that the particle dynamics is compatible with the rotation 
group SO (3) since the velocity field is constrained to be an element of the Lie 
algebra so(3) of the rotation group, i.e., Vi = VLijXj, where the 3x3 matrix 
VLij = eijkOJk is an antisymmetric matrix. Recall that the set of all 3 x 3 real 
antisymmetric matrices forms the Lie algebra of the rotation group. 



The motion group for the quantum Bohr-Mottelson collective model M 



and the classical Riemann ellipsoidal model |11| is the general linear group 



GL(3,R). The general linear group constraint means that the velocity field of 
the particle system is a linear function of the Cartesian position coordinates, 
i.e., Vi = '^ijXj, where S is a constant 3x3 real matrix. Thus S is an element 
of the Lie algebra gl(3,R), the set of all 3 x 3 real matrices. Rigid rotation 
requires that H = i7 be antisymmetric, while irrotational flow (V x -u = 0) 
demands that S be symmetric. Divergence- free vector fields {V ■ v = 0), 
corresponding to incompressible flow, imply that S is traceless, i.e., S is an 
element of the Lie algebra sl(3,R) of the special linear group. 

A serious limitation imposed on these familiar geometrical models is that 
the velocity fields are always linear. This strong condition is violated unam- 
biguously in some circumstances, e.g., fissioning isotopes with a neck or other 
exotic shapes, and possibly broken strongly for other nuclear states, even at 
low energy. The article has two interrelated objectives: first, to determine 
nonlinear vector fields on Euclidean space that, for example, may model the 
dynamics of a fissioning isotope, and, second, to construct collective models 
compatible with such nonlinear motion groups. 

The first problem is a surfeit of riches: there is an K-infinity of nonlinear 
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vector fields from wliicli to cfioose tlie model's collective dynamics. Indeed 
every triple of smooth functions on Euclidean space defines a different vector 
field. Among all possible vector fields, a subset V of collective motion vector 
fields will be chosen. To be relevant to a dynamical system of interacting 
particles, the subset must assimilate the tangents to the particle trajectories; 
to be useful, these collective vector fields should not include much more. 

The collective motion vector fields satisfy three general properties. First, 
V contains the angular momentum algebra so(3). Because of the SO (3) in- 
variance of the nuclear Hamiltonian, rotated states are degenerate in energy. 
Thus collective rotational motion is expected as a rule. An exception to this 
typical situation occurs when rotational symmetry is broken, say by an ex- 
ternal magnetic field. In this special case the set of collective vector fields V 
may contain only so (2). 

Second, the collective velocity fields must form a finite-dimensional set. 
By definition, "collective" particle motion means that the particle velocity 
vectors are constrained. Indeed collective motion requires that the particles 
move together in an orchestrated manner, and, therefore, the physically ad- 
missible particle velocity vectors are restricted to a finite-dimensional set. 
If the set V were infinite dimensional, then no collective physics would be 
infused into the model's ansatz. At best, the complexity of the model would 
be similar or identical to that of the original many-body problem. As a 
practical matter, the dimension of the set of collective vector fields should 
be much less than the dimension of the A-particle Euclidean configuration 
space, dimV << 3 A. 

Finally, V is a Lie algebra, i.e., the set of collective vector fields is closed 
under commutation. One reason is that the family of collective motions 
is closed under composition. If X is a collective vector field, let exp(tX)x 
denote its integral curve through the point x. If X and Y are two collective 
vector fields, then an admissible collective motion through x is given by the 
composition of collective modes 

7(i) = exp{ViX) ex.p{ViY) ex.p{-ViX) ex.p{-ViY)x. (1) 

Since the tangent to the curve 7 is the vector field [X, Y] , the commutator 
of two collective vector fields must also be a collective vector field. 

Another reason is that this algebraic property is necessary to formulate 
a quantum model of collective motion corresponding to the geometric vector 
fields in V. The natural construction of the quantum theory from the Lie 
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algebra of collective vector fields is given by a series of well-defined mathe- 
matical operations; the archetypes for the construction are the microscopic 
theories corresponding to the adiabatic rotational and Bohr-Mottelson mod- 
els Ijl^, 0, |14|, |15|, |1^. First, the vector fields on Euclidean space lift to 
vector fields on the A-particle configuration space R^"^ by simply summing 
over the particle index to create one-body operators. We seldom distinguish 
between V, considered as vector fields on R^, and the isomorphic algebra of 
vector fields on R^"^. The sum over all particles insures that both the Pauli 
principle is obeyed and that the generated motion is collective. 

Second, according to Frobenius's theorem, the collective vector fields on 
R'^'^ may be integrated to surfaces in the A-particle configuration space |17 . 
The vector fields in V are tangent everywhere to these integrated surfaces. 
Each surface is an orbit of a connected Lie group G whose Lie algebra is V. 
Third, each orbit is a homogeneous space diffeomorphic to the coset space 
G/H, where H is the isotropy subgroup of the orbit surface. Excluding a set 
of measure zero, the isotropy subgroup is the same for all orbits; the excep- 
tions are orbits of smaller dimension than the generic orbit. The common 
dimension of the generic orbit surface is the difference between the dimen- 
sions of G and its subgroup H, dimG/H = dim G — dim if. If the dimension 
of the y4-particle configuration space is greater than the dimension of V, 
3 A > dimV, then the generic isotropy subgroup H is expected typically to 
be discrete. 

Fourth, a transversal submanifold Af of R^"^ is chosen that intersects 
each generic orbit orthogonally exactly once. The transversal manifold is 
not unique. After removing a set of measure zero (the lower dimensional G- 
orbits), the A-particle configuration space is diffeomorphic to the Cartesian 
product of G/H and the transversal manifold A/". Since sets of measure zero 
are irrelevant to the Hilbert space of square-integrable measurable functions 
on R^^, this Hilbert space is isomorphic to the tensor product, 

C\R^^) = C\G/H) ® C^{U). (2) 

The collective vector fields of V act only on the first component C^{G/H) of 
the tensor product; this space is the Hilbert space of collective wave functions. 
Intrinsic wave functions are vectors in the space C'^{Af). Bands of collective 
states are given by monomials (PtKim ® X; where x is a fixed intrinsic state in 
the space of the transversal submanifold and (prKiM are bands of collective 
states from the space of the homogeneous manifold G/H. 
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A familiar example illustrates this somewhat abstract construction: let 
V be the angular momentum algebra so(3) and suppose there is only one 
particle A = 1. In this case the collective vector fields integrate to the 
surfaces of spheres S2', the angular momentum vectors are tangent to the 
spheres. Each sphere is an orbit of the rotation group SO (3) and the sphere 
is diffeomorphic to the coset space S2 = S0{3)/ S0{2). An exception is 
the singular orbit consisting of a single point, the origin. An orthogonal 
transversal is any straight line through the origin, e.g., the positive 2;-axis. 
Thus the Hilbert space of square-integrable wave functions is isomorphic to 
the tensor product 

£2(R3) = C\SO{3)/SO{2)) ® C\R+). (3) 

A basis of "collective" wave functions consists of the spherical harmonics; 
intrinsic states are just radial wave functions. 

For a successful decomposition into collective and intrinsic wave functions, 
it is absolutely essential that Frobenius's theorem may be applied. But the 
necessary and sufficient condition for the existence of an integrating surface 
is that the set of collective vector fields close under commutation. Hence 
there is no fiexibility here: V must be a Lie algebra. 

The mathematical problem of classifying all nonlinear vector fields on Eu- 
clidean space that close to form a finite dimensional Lie algebra was studied 



first by Sophus Lie |]TB[. In one dimension there are only three possibilities: 
dilations (ax), affine maps {ax + b), and linear fractional transformations 
{{ax + b)/{cx + d)). In two dimensions there is already a large number of 
possible nonisomorphic Lie algebras of vector fields. In three dimensions, be- 
cause the number is so large, there is no known enumeration of all solutions. 
Lie's classification, however, is at once too fine and too coarse. It is too fine 
because we restrict to Lie algebras of vector fields that contain the angular 
momentum. It is too coarse because Lie does not distinguish among vector 
fields that differ by changes of variables: these are isomorphic Lie algebras 
of vector fields. But physics does distinguish among velocity vectors that 
differ by point transformations since the particle trajectories are different. 
Indeed physics distinguishes among Hamiltonians that differ by canonical 
point transformations. Thus Lie's classic work is not directly relevant to our 
physics problem. 

As an initial investigation of nonlinear collective nuclear motion, the 
smallest simple Lie algebra of smooth vector fields on Euclidean space that 
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contains the angular momentum algebra will be determined. This minimal 
extension of so(3) is isomorphic to the algebra of the special linear group, 
V ~ s/(3,R) D so{3). It will be shown in section two that the set of all 
sl(3,R) nonlinear vector fields is indexed by the set of smooth real-valued 
functions in one real variable. If each nonlinear vector field is also divergence 
free (because the motion is incompressible), then the set of incompressible 
nonlinear sl(3,R) vector fields is indexed by a single real number A. 

Among the solutions is the usual linear representation of sl(3,R) corre- 
sponding to A = 0. The linear vector field is defined by the Lie derivatives 
of the natural linear action of the group SL(3,R) on three dimensional Eu- 
clidean space. In section three it is proven that all incompressible nonlinear 
sl(3,R) vector fields are given by certain nonlinear point transformations ap- 
plied to the linear representation. The spaces of square integrable functions 
in the nonlinear geometries (A ^ 0) are shown to be isometric to the usual 
Hilbert space of square integrable wave functions (A = 0). 

In section four the algebra of nonlinear sl(3,R) vector fields is extended to 
the Lie algebra of gl(3,R). A family of surfaces is discovered that is invariant 
with respect to the nonlinear GL(3,R) motion group. These surfaces play 
the same role in the physical interpretation of the nonlinear theory as the 
inertia ellipsoids do in the linear case. The nonlinear surfaces include shapes 
with a neck (A > 0) and also configurations with a bubble ( A < 0) . 

In section five the construction of nonlinear collective models is achieved. 
A closed fifteen dimensional Lie algebra, isomorphic to gcm(3), is presented 
that contains both the gl(3,R) algebra of the nonlinear motion group plus 
certain nonlinear quadrupole and monopole tensors. When A = the non- 
hnear gcm(3) algebra simphfies to the usual linear gcm(3) algebra of the Rie- 
mann ellipsoidal model. All nonlinear collective models compatible with the 
nonlinear GL(3,R) motion group are given by the irreducible unitary repre- 
sentations of the GCM(3) Lie group. These representations are determined 
completely by the inducing construction because GCM(3) is a semidirect 
product with an abelian normal subgroup. 

Because of the relevance to the nonlinear theory of the natural linear 
representation of the general linear group on three dimensional Euclidean 
space, a synopsis of its properties is provided here to close the introduction. 
The general linear group contains the familiar rotation group SO (3) as a 
subgroup. In addition, GL(3,R) contains the subgroup of invertible diagonal 
matrices that transform space by expanding and compressing the Euclidean 
X — y — z axes lengths. Yet another subgroup is the nilpotent group of up- 
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per triangular matrices with ones running down the diagonal; these group 
elements are the shear transformations. The special linear group SL(3,R) 
of unimodular matrices is a subgroup of the general linear group; SL(3,R) 
consists of the incompressible linear transformations [0, An important 
property of the general linear group is that its action transforms an ellip- 
soid into another (deformed and rotated) ellipsoid. This property is at the 
foundation of the Riemann ellipsoidal model ^ . 

The unitary representation U of the general linear group on the Hilbert 
space of A-nucleon wave functions is provided by 

{U{g)^){n,r2, . . . ,r4) = {det gY^^^ig-' ■ n, g~' ■ r^, ...,g-'- ^a), (4) 

where g G GL(3,R). The quantum mechanical hermitian observable that 
corresponds to each one-parameter geometrical group action on Euclidean 
space is defined by the Lie derivative of the representation U. The angular 
momentum operators, Lij = XiPj — XjPi, summed over the particle index, 
span the Lie algebra of the rotation group. The breathing mode oscillation 
operator 5* is the Lie derivative of the GL(1,R) unitary representation of 
uniform scaling transformations 

S = f-p- {3/2)ih, (5) 

summed over the A-particles. For the transformation g = diag(e'^, e"^, e*^) 
of Euclidean space, the corresponding unitary operator is given by U{g) = 
exp{—ieS/h). 

The SL(3,R) subgroup generates the angular momentum and a second- 
rank tensor operator T, 

Tij = XiPj + XjPi - (2/3)f ■ p Sij, (6) 

where 1 < "i, j < 3 label the Cartesian axes x, y, z and the expression must be 
summed over the A-particles. To attain a geometrical understanding, note 
that the unimodular transformation 

g = di8.g{e-%e-%e^^) (7) 

corresponds to the unitary operator U{g) = exp(— 3ieT33/2/i). There are 
similar expressions for incompressible scaling in the x and y directions. The 
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off-diagonal components of the symmetric tensor Tj^ are associated with the 
hyperbola mappings 



for which U{g) — exp{—ieTi2/h). 

2 Nonlinear Vector Fields on 

In this section all smooth vector fields on three dimensional Euclidean space 
are determined that (1) obey the sl(3,R) commutation relations and (2) 
contain the standard angular momentum so (3) subalgebra. The results are 
summarized in the theorem at the end of this section. The linear vector field 
representation of sl(3,R), as defined by the Tj^'s above, is one acceptable 
solution. Is it the unique solution to conditions (1) and (2), or are there 
other nonlinear solutions? 

Were it not for the so (3) subalgebra limitation, the exhaustive compu- 
tation of all sl(3,R) vector fields would be a formidable task with many 
solutions. When the so(3) constraint is imposed, the set of all smooth non- 
linear sl(3,R) vector fields will be proven to be indexed by a single smooth 
function A{r) of the radial coordinate r. The linear solution is just one of 
the possibilities, A{r) = r. If the sl(3,R) vector fields are restricted further 
to have zero divergence, then it will be shown that A{r) = r + A/r^ where A 
is a real constant. 

Because the so(3) subalgebra is fixed, spherical coordinates and bases 
are most convenient for the calculations. The eight dimensional simple Lie 
algebra sl(3,R) is spanned by the vector angular momentum L and a second 
rank tensor T. The angular momenta span so (3), the maximal compact 
subalgebra of the unimodular algebra. The angular momentum vector fields 
are given in spherical coordinates (r, ^,0) by 





.d_ 



= e-^[±- + .cot.-]. 

These vector fields obey the so(3) commutation relations 

[Lo,L±] = ±L± 



(9) 
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(10) 



In the spherical basis, the commutation relations require that the com- 
ponents form a second rank tensor operator, 

[Lo,T^] = /iT^ 



[L±,T,] = v/6-^(/i±l) Vi, (11) 

for /i = — 2, . . . , +2, and that the set closes to form the simple Lie algebra 
sl(3,R), 

[T^,T,] = -2v^(2/^2i/|l^ + z/)L^+„ (12) 

where L±i = ^{1/ \^)L±. The spherical basis is normalized to the Cartesian 
basis by 

T±2 = ^(Tn-T22)±^Ti2. (13) 

In particular, the spherical component T+2 of the linear vector fields of 
Eq. (^, is given by 

= -.e-^[rsin^e|- + ^sin2^^ + ^^]. (14) 

To find all possible sl(3,R) vector fields, we start with a completely general 
form for T+2, 

o o o 

T+2 = [a(r, OA)g^ + b{r, OA)g^ + c(r, 9, 0)—], (15) 

and then impose the commutation relations upon it. The coefficients a, b, c 
of the vector field are assumed to be smooth functions of the spherical co- 
ordinates in an open subset of Euclidean space that is contained within the 
domain of the spherical coordinate chart, viz., r>0, 0<^<7r, O<0< 27r, 
or minus the 2;-axis. In addition each coefficient is periodic in (p with 
period 2n. 

First, [Lo,T+2] = 2T_|.2 determines completely the dependence of the 
coefficients, 

T^2 = -^ e^'Mr. 9)^ + Kr, 0)^ + c(r, 9)^]. (16) 
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Second, the highest weight condition, [L+, = 0, requires that three inde- 
pendent equations be satisfied, viz., the vanishing coefficients of the partial 
derivatives in the r, 9, and directions, 

— {r, 6) -2 cote a{r,e) = (17) 

db 

— {r, 9) ~ 2 cote b{r, 9) -ic{r, 9) = (18) 

CjC 

— (r, 9) - cot 9 c(r, 9) + i csc^ 9b{r,9) = 0. (19) 

Eq. (|l^) can be solved immediately to give 

a(r, 9) = A{r) sin^ 9. (20) 

If 6(r, 9) is expressed in terms of c(r, 9) using Eq. (^) and substituted into 
Eq. (|l^), then the resulting second order differential equation may be solved 
for c{r,9). The general solutions are 

b{r,9) = i sin 9 [a [r) — j3{r) cos 9] 

c{r,9) = -a{r) cos 9 + p{r) (21) 

where A{r), a{r), and j3{r) are arbitrary smooth functions of r. In addi- 
tion, repeated application of the lowering operator L_ to the highest weight 
determines the lowest weight component of the tensor T, 

T±2 = -ie^^''''[A(r)sm'^9— ±i sin 9 (a(r)Tf3 (r) COS 9) — 

or o9 

±{(3{r)Ta{r)cos9)^]. (22) 

Because of the Jacobi identity and the fact that T is a tensor operator, it is 
sufficient for the one commutation relation, 

[T+2,T_2] = -4Lo, (23) 

to be satisfied to guarantee that all the commutation relations among the 
components of T, Eq. (|12D, are correct. Eq. (|23| ) yields three equations, 

= A{r)a{r) (24) 

= A(r)^-3ia(r)p(r) (25) 
dr 

4i = -4i(3{ry + 4ia{rfcosH + 2(^A{r)^-ta{r)^^sin^9. (26) 
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Since Eq. (^Bj) must be satisfied for all < < tt, the coefficients of cos^ 6 
and sin^ 6 are zero while the constant term on the right hand side must equal 
Ai. The solution is a(r) = and /3(r)^ = —1 or /9(r) = ±i. With these 
choices for a and /3, Eqs.(0,^) are satisfied too for any function A{r). If 
/3 = +z, then assembling our results, 



T±2 = -ie^'^'^[A{r) sin' 0^ + - sin 26^ ± i 



d 



d 



d 



dr 



de 



(27) 



for some smooth function A{r) of the radial coordinate. If A{r) = r, then 
the linear representation of sl(3,R) is recovered. 

From Eq. (^) it is evident that, if is a solution to the sl(3,R) commu- 
tation relations, then a multiple cT^ is also a solution if and only if c' = 1 
or c = ±1. Note that (3 = ±i corresponds to c = ±1. For simplicity, the 
ambiguous phase was chosen in Eq. ( ]27|) to match the conventional phase of 
the linear vector fields, Eq. (0). 

Thus the set of smooth sl(3,R) vector fields on that contain the usual 
so (3) subalgebra of angular momentum vector fields is indexed by a smooth 
function A{r). We now see how strict the angular momentum constraint 
really is. All the angular dependence of the T^'s has been determined and 
the only freedom left is in the radial component. 

The radial component may be constrained further if the motion is incom- 
pressible. In this case, the sl(3,R) vector fields must be divergence free. But 
the divergence for the highest and lowest weight components is 



div T-t2 = —i e 



The divergence vanishes if 



±2i4 



l_d_ 



r'y4(r)) 



(28) 



A{r) 



r + 



A 



(29) 



where A is a real number with the dimensions of volume. Summarizing the 
results, we have proven the following: 



Theorem. For every smooth function A{r) of the radial coordinate r, 
there is a representation of the simple Lie algebra sl(3,R) by smooth vector 
fields in an open domain of R^ such that (1) the so(3) rotation subalgebra 
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is represented by the standard angular momentum operators and (2) the 
components of the rank two spherical tensor are given by 



o o o 

T±i(A) = -ie^^^[TAir)sm2e—Tcos29—-icot9—] 
To{A) = Zl[A(r)(l + 3cos2^)— -3sin20— ]. 

If the representation is restricted to divergence free vector fields, then A{r) = 
r + h./r'^, where A is a real number. Note that A = is the usual linear 
representation. 

As a corollary, there is no similar representation of the compact Lie alge- 
bra su(3). The reason is that, since the complexifications of sl(3,R) and su(3) 
are isomorphic, the theorem determines every solution to the su(3) algebra 
too, viz., iT^. This su(3) tensor corresponds to a "vector field" with pure 
imaginary components. 



3 Origin of nonlinear sl(3,R) vector fields 

In this section every divergence-free so(3)-preserving nonlinear representation 
of sl(3,R) is shown to be given by a change of variables applied to the linear 
representation. The proof is achieved using the concept of a manifold. 



Definition. A set M is a smooth manifold if: ||T9| 



1. M is a topological space with a family of open sets Mq, that cover it, 

2. Associated with each such open set is a one-to-one bicontinuous map- 
ping ifja from Mq onto an open set in R". Each pair (M„,'?/^q,) is 
called a chart and the set of all charts is called an atlas for the manifold 
M. 

3. Given two different charts, (Mq,, V'o) and (Ma,'?/'/^), the map ipp o 
from the open subset ipa^Ma fl Mg) of R" to the open subset iljpi^Ma n 
Mp) of R" is smooth, i.e., infinitely differentiable. 
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It is important to maintain the distinction between a manifold and the charts 
that provide its coordinatization. For the Euchdean space M = R3 this 
distinction is often ignored because the usual atlas, Cartesian coordinates, 
consists of only a single chart (Mxyzjipxyz), where M^yz = Oxyz = R-^, and 
ipxyz : m I— > (x, y, z) is the identity mapping for every point m e R^. Another 
common chart for three dimensional space is given by the spherical coordi- 
nates. But spherical coordinates do not supply a global chart for because 
the homeomorphism iprecf, : m {r, 6, 0) is not one-to-one everywhere. Along 
the 2;-axis, 9 = and 9 = ti, all values of the longitudinal angle 4> correspond 
to the same point of the manifold M = R^. The situation is even worse at 
the origin, r = 0, where all angular coordinates 9 and (j) correspond to the 
same point of the manifold. In addition, since the longitudinal angle range 
[0,27r] is a closed set, at least two overlapping charts are required to cover 
R^ minus the z-axis. For example, one chart (M^^o, 'V're^) excludes = 0, 
while the other (M^^^, ^^ec/i) excludes = tt. Although the chart mapping 
ipre(i> is the same for both regions, the domains and images of these charts are 
different: 

O^^o = ^re^iM^^o) = {(r, ^, 0)|r > 0, < e < TT, < < 27r} 
O^^^ = ^,9<^(M<^^^) = {(r,^,0)|r>O,O<^<7r,-7r<0<7r}. (30) 

Evidently similar charts (M^^j^^, ■j/'re^) may be defined by excluding an arbi- 
trary longitudinal angle 0o from the range O^^ff,^. Let {Mro^, iprdrf)) denote 
any such spherical coordinate chart mapping onto the open set ipre<i>{Mrd(j,) = 
OrOtji- The third condition in the definition of the manifold is satisfied by the 
smooth function ^l)xyz ° "07^0 ■ ^rO^) ^ tpxyziMrOcf)), where x = rcos0sin^, 
y = r sin sin 9, z = r cos 9. 

As will be proven shortly, the incompressible nonlinear vector fields are 
described naturally by two charts {M^,ip^) and(M_,'0_) corresponding to 
A > and A < 0, respectively. The open set M+ = is identical to 

the domain of a spherical coordinate chart, but the range = ■0+(M+) 
excludes a sphere of radius R = \/A: 

^+ ° ^r"i : ^re<p — ^ 0+ 

^r^OA) ^ (r',^',0') = (v/^^,^,0), (31) 

where (9+ = OyQ^ fl Sr and Sr denotes the complement of the closed sphere 
of radius R. If A < 0, the domain of the chart M_ = M^g^ fl S'r is a region 
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excluding the sphere of radius R — but now the range is Ore4> — 

ijj-{M^) (see Figure 1): 



(^,^,0) ^ (r',^',0') = (Vr^,^,0). (32) 

The next step is to demonstrate that the nonhnear vector fields of §2 can 
be constructed in a simple way via a change of chart applied to the linear 
vector field representation. Suppose {Ma,ipa) is any chart of the manifold 
and Oa — ipa{Ma) is an open subset of R^. The Lie algebra sl(3,R) is 
represented by linear vector fields on the open set Oa- In particular, if 
spherical coordinates (r', 9\ (j)') are used for the open set Oa, then the highest 
weight component of the rank two sl(3,R) vector field is given explicitly by 

= -ze-^'[r'sin^^'A + ^ sin 2^' A + z A]. (33) 

A vector field defined on is mapped back to in the Euclidean 

space by ip'^^ and then expressed in terms of the usual spherical coordinates 
by V^r6»(}!) ° V^a ^- Hcncc, the Euclidean space vector field, corresponding to 
on is 

(IT 1 r) ff 

= -e-[(r+^)sin^.|+l,.,2.| + .|] 

= T+2(A), (34) 

where {ipre(i> ° '^+^)* denotes the Jacobian of the transformation iprecp ° i^a^- 
Indeed every component of the nonlinear sl(3,R) representation for A > 
defined on the open set iV'/_|_ = M^g^ is given by the coordinate change 

T,{K)^{^,e4>o^l%Tl. (35) 

If A < 0, a similar argument demonstrates that the nonlinear representation 
on the open subset M_ = M^g^ fl S of R^ is given by the above equation 

when ^_|_ is replace by ip-. Thus, for A < 0, the nonlinear representation is 
defined in a region of the Euclidean space exterior to the sphere of radius 
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Note that, for any A, the so(3) subalgebra is preserved by the mappings 

L^ = {i)re,po^^^)^L'^. (36) 

It is now clear why the nonlinear vector fields are a representation of the 
sl(3,R) algebra: The Jacobian mapping {ilJre<p°'ipa^)* is a Lie homomorphism 
of the algebra of linear vector fields defined on 0^ into the algebra of, in 
general nonlinear, vector fields on OrOcj, [0 • 

The nonlinear vector fields may not be complete. The linear vector fields 
are certainly complete on R^, i.e., they can be integrated to one-parameter 
subgroups of the Lie group SL(3,R), e.g., Eqs.( |^, || ). If a one-parameter 
subgroup is contained entirely within the region (9^, then the nonlinear vector 
field is complete because the subgroup is mapped entirely within by ip~^ ■ 
Otherwise the vector field may not be complete. For A < 0, the nonlinear 
vector fields on M_ are complete since C_ = Ore<p- Hence, SL(3,R) acts on 
M_ as a Lie transformation group. But for positive A the region 0+ excludes 
the sphere of radius v^, and integral curves of the linear representation will 
typically run into this sphere. In this case, the integral curves of the nonlinear 
sl(3,R) vector fields form one-parameter semigroups instead of one-parameter 
groups. 

Thus the connection between linear and nonlinear vector field represen- 
tations of sl(3,R) is established. To complete the picture, wave functions in 
different coordinate charts must be related to each other by reformulating 
the chart mappings as isometries between different Hilbert spaces. 

If M is a manifold and o; is a volume element on it, then define the Hilbert 
space £^(M) to be the set of all measurable complex- valued square-integrable 
wave functions on M: 

C\M) = |/ : M ^ C I / is measurable and^^ |/pcu < ooj . (37) 

A wave function / is expressed in a chart {Ma,ipa) as a complex- valued 
function fa = f ° i^a^ '■ Oa C. Given two overlapping charts {Ma,ipa) 
and [M/Sjipp), the coordinate expressions for the wave function / are related 
by a change of variables, f/3 = fa° ii^a ° 0/?^)- The expressions for the 
volume element in the two charts are related by the determinant Jajj of the 
Jacobian of the coordinate transformation %jja o 0^^,ci;q = Jafsujp. If a single 
coordinate chart covers the entire manifold with the possible exception of a 
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set of measure zero, then the Hilbert space £^(M) is isomorphic to the space 

C^{Oa) = |/a : Oa — > C | is measurable and j \fa\^uJa < ooj , (38) 

where uJa denotes the volume element uj as expressed in terms of the chart's 
coordinates. More generally, wave functions and integration on the manifold 



M are pieced together from different charts using a partition of unity p9 |. 

For example, the Hilbert space C'^{R?) of square- integrable wave func- 
tions on the three dimensional Euclidean manifold is isomorphic to the space 
of complex-valued measurable wave functions of the Cartesian coordinates 
{x,y,z) that are square-integrable relative to the volume element dxdydz. 
Since the 2;-axis is a set of measure zero, the space of complex-valued measur- 
able functions of the spherical coordinates (r, 6', 0) that are square-integrable 
on the region Or^t^ relative to r'^ sin 6 dr d6 d(j) is also isomorphic to £^(R^). 

Since M+ = MrOcf,, the space is also isomorphic to £^(R'^), 



(39) 



£^(C+) = |/+ : (9+ ^ C I /+ is measurable and J |/+|^ti;+ < oo| 

The volume element in the chart coordinates (r', 6', (J)') is computed from 
the Jacobian of the transformation from the spherical Ore^ chart coordinates 

a;_|_ = r'^dr sin 9d9d(j) 
dv 

= r'^— —dr' sin 9 dOdcj) 
dr' 

= {r' f sine' dr'dd'd<p'. (40) 

Since u;+ is the same as the spherical coordinate volume element, the nonlin- 
ear vector fields T^(A) for A > form a hermitian rank two tensor operator. 

Each Hilbert space £^((9+) for positive A is therefore canonically isomor- 
phic to the Hilbert space £^(R'^) = C?{Ore^). The explicit isometry is given 
by 

UkJ = f o (^ip^Q^ o ^^1) 

(f/A/)(r',^,0) = f{^?^^,e,(l>). (41) 
The isometrics obey the additive property 

Ua.Ua, = Ua,+a„ (42) 
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and the inverse mapping is f/_A : C'^{OrB,j,)- 

The infinitesimal generator for = exp(— iAG) is given by 

d 



= -^ir.%ir'A^). (43) 

G is hermitian because U\ is an isometry. In terms of the hermitian radial 
momentum operator, = —ir'^{d/dr)r, the generator is 

G = ^p^. (44) 

Sr r 

The nonlinear sl(3,R) operators may be expressed as a sum of the linear 
representation plus the product of the rank-two spherical harmonics times 
the generator G, 

T,(A) = T,(0) + ^AY,G. (45) 

The situation for negative A is more subtle because the chart (M_, ■?/'_) 
omits the sphere of radius R = |A|^/^ from the Euclidean manifold. In gen- 
eral, a second chart is required to cover the interior of the sphere. However, 
in the special case of bubble nuclei, the single-particle wave functions must 
vanish inside a sphere of radius R. The subspace of bubble nuclei wave 
functions is isomorphic to 

/:2(C»_) = I C»_ ^ C|/_ is measurable and ^ < ooj . (46) 

The volume element U- is again the same as the spherical coordinate volume 
element, and, hence, the nonlinear vector fields T^(A) for A < also form a 
hermitian rank two tensor operator. The isometry U\ for negative A is given 
by Eq. (^T]), but the range is now £^((9_). 

The sl(3,R) results may be extended easily to nonlinear gl(3,R) vector 
fields. The hermitian monopole operator of gl(3,R) is given in a coordinate 
chart {M±,ip±) by 

S' = -zr'^ - -I. (47) 
dr' 2 ^ ^ 
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On the Euclidean manifold, the corresponding nonlinear vector field is 

. ,dr d 3 . 

dr' dr 2 

= 5(0) + SAG. (48) 



4 Deformed Droplets 

Consider an ellipsoid in three-dimensional Euclidean space. Its boundary, 
defined uniquely by a positive-definite symmetric 3x3 real matrix consists 
of the points that satisfy the quadratic relation 

(49) 

where r denotes the Cartesian position vector. If g is a diagonal matix with 
entries a^, 6^, and c^, then the ellipsoid's principal axes are aligned with the 
Cartesian axes and its semi-axis lengths are a, 6, and c. 

With respect to the usual action of the general linear group, r — > gf for 
g G GL(3,R), an ellipsoid is transformed into another deformed and rotated 
one, q — > g-q-^g. This action transforms the sphere of unit radius {q = I) into 
the ellipsoid corresponding to the positive-definite matrix q — g-^g. But any 
positive-definite matrix may be expressed in the form g ■ *g. Therefore, the 
general linear group acts transitively upon the family of ellipsoidal surfaces. 

With respect to the nonlinear actions considered in this article, GL(3,R) 
also acts transitively upon a family of surfaces, but these are not ellipsoidal 
boundaries anymore. As shown in §3, the nonlinear group action on open 
subsets M± of the Euclidean manifold is constructed by applying the nonlin- 
ear coordinate transformation ip:^^ to the linear representation of GL(3,R) in 
the coordinate chart 0±. Hence the invariant family of nonlinear surfaces in 
the Euclidean manifold is given by transforming ellipsoids in the chart space 
0± back to the manifold via t/jg^. If Cartesian coordinates are used for the 
Euclidean manifold, then each nonlinear surface in the invariant family is 
given by the solutions to 

/ A \ 

1 = tf' g-i f ' = f 1 + —J *f g-^ f, (50) 

where f ' and r denote Cartesian coordinates for 0± and the Euclidean man- 
ifold, respectively. 
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If A = the figures are ellipsoids. If A is negative, then the ellipsoidal 
solutions to Eq. (0) are contained entirely within the chart space C_ = Orecj,- 
Its pre-image in M_ is the surface of a deformed droplet excluding a spherical 
hole of radius \/|A|, see Figure 2. The half-lengths of the droplet's principal 



axes are i^a^ + |A|, -^6^ + |A|, and y + |A|. The volume of the droplet in 
Euclidean space, excluding the spherical bubble of volume 47r|A|/3, equals 

the volume Anabc/S of the ellipsoid in the chart space. For ^/|A[ 3> a, b, c, the 
nonlinear droplet is a thin shell surrounding a spherical bubble. The width 
of this shell is not uniform; along a principal axis, say the x-axis, the width of 
the thin shell is approximately (a^/|A|)^/^a/3. In addition to exotic bubble 
nuclei, the negative-A figures may apply to the invariant core or two-fluid 
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If A is positive, then the ellipsoidal solutions to Eq. may not be 
contained entirely within the chart space (9+ = Ore<i, fl Sa^. The pre-image 
of the portion of the ellipsoid contained in the interior of the sphere of radius 
is the empty set. The intersection of the ellipsoid's surface with the 
sphere's surface is mapped by ip^^ to the Euclidean space origin. The pre- 
image of the ellipsoidal boundary exterior to the sphere is the surface of a 
deformed droplet in the Euclidean manifold. In Figures 3 and 4, several 
positive- A nonlinear surfaces are drawn for axially symmetric (a=c=l, b=2) 
and triaxial (a=l, b=3, c=2) shapes, respectively. For the axially symmetric 
solutions, the nonlinear figures are surfaces of revolution that develop a neck 
as A increases. At v^A = a = c = 1, the neck is pinched down to a single point 
because the sphere in the chart space touches the two short principal axes 
of the prolate ellipsoid. For A > 1, the droplet's pinched neck is elongated, 
until, when = 6 = 2 (the length of the long axis of the prolate ellipsoid), 
the figure collapses to a point. Beyond that there is no figure in Euclidean 
space. 

In Figure 4, the surfaces of several triaxial deformed droplets are drawn. 
As advances through the half-lengths of the principal axes a = 1, c = 2 
and b = 3, the pinched neck first is squeezed in the x-direction, then collapses 
to a point, and finally the figure disappears altogether. The volume of the 
Euclidean space droplet is the difference between the volume Airabc/S of the 
ellipsoid in the chart space and the portion of the volume of the sphere of 
radius contained in the interior of the ellipsoid. In particular, if \/7i 
is smaller than any of the half-lengths, a, b, and c, then the volume of the 
deformed droplet equals 47r(a6c — A)/3. Some of the shapes in Figures 3 and 
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4 are similar to the equilibrium surfaces of a rotating charged droplet near 
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To determine the mechanics of a continuum fluid or finite particle sys- 
tem corresponding to a deformed droplet, an observable must be identified 
that characterizes these nonlinear geometrical shapes. Consider first a fluid 
contained within an ellipsoidal boundary with uniform density po- Ifs axes 
lengths and orientation are determined by its quadrupole plus monopole Q 
tensor 

Qij = J Po(r)xiXjd^x. (51) 

The inertia tensor for a rigidly rotating system is related to the Q tensor by 
lij = TT{Q)6ij — Qij. If the principal axes are aligned with the Cartesian 
axes, the Q tensor is diagonal with entries (M/5)a^, (M/5)6^, and (M/5)c^, 
where M denotes the ellipsoid's mass. The orthogonal transformation that 
diagonalizes the Q tensor is the rotation that aligns the Cartesian axes with 
the ellipsoid's principal axes. Hence the Q tensor characterizes an ellipsoid's 
size and orientation. 

If A is negative, the surface of the deformed fluid droplet in the Euclidean 
manifold is mapped by ip- to an ellipsoidal surface in the chart space 0-. 
Suppose the density of the droplet in the Euclidean space is given by p = 
Po o "ip- on M_ and equals zero in the interior of the sphere of radius '^flM- 
Then the Q{A) tensor is defined by 



Qij{A) = I po{r)x'ix'jd x' 

A 



j p{f){l + -f'^XiXjSx, (52) 



where the first integral is over the chart domain = OrOij, and the second 
is over the Euclidean manifold excluding the sphere of radius ^|A|. If the 
ellipsoidal surface is defined by the matrix g, then Q(A) = (M/5)g. 

If A is positive, the deformed droplet's surface is mapped by ip+ to, in 
general, a portion of an ellipsoidal surface in O^. The excluded portion is 
the intersection of the ellipsoid with the sphere of radius v^. Suppose the 
density of the deformed droplet in Euclidean space is p = poo ■?/;_,_. Then Q{A) 
is defined again by Eq. (p2D where the first integral is over the domain Oj^ 
that excludes the sphere and the second is over the entire Euclidean space. 
But, the correspondence between Q{A) and the matrix q that defines the 
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boundary of the ellipsoid is not so simple. If is smaller than any of the 
half-lengths, a, b, and c, then the first integral of Eq. (0) is evaluated to be 

,53) 

For a discrete classical system of A identical particles the Q tensor sim- 
plifies to a sum over the particles 

A ^ 

Q.j(A) =E (1 + -f^'^mXnj, (54) 
n=l 

where the mass density is replaced here by the number density. To recover 
the mass density expression, merely multiply Eq. (^^ by the common mass 
of the particles. 

For a quantum system of A identical particles described by a wave func- 
tion the expectation of the one-body operator Qij(A) of Eq. (^41) equals 
the integral Eq. (|5^) where the one-body density is 

p{f) = A j \^{f,r2,...fA)\'^d^r2---(frA. (55) 

The monopole and quadrupole moments of the deformed droplet are the 
quantum expectations of the one-body operators 

Qo(A) = tirl + ^f" 

n=l 

Q2mW = E(r3+A)2/3F2^(^„,0„), (56) 



n=l 



where Y2m.{0,4>) denotes the mth component of the rank two spherical har- 
monic. If r„ ^ \f\M^ then the monopole, quadruple, and Q tensors ap- 
proximate their usual values. The Bohr-Mottelson parameters (/3, 7) for the 
deformed droplet may be defined in terms of the expectations of the quadratic 
and cubic scalars |2^, |2^, ^ 



< [Q2(A) X g2(A)]° > oc 
< [Q2(A) X Q2{A) X g2(A)]° > oc /^^ COS 37 (57) 
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The square of the volume of the deformed droplet is related to the expec- 
tation of the determinant of Qij{K). For an ellipsoid (A = 0), the square of 
the volume is proportional to the determinant of the Q tensor 

^^-Q)'(f)'<detQ>=0. (58) 

If A < 0, then the left hand side of Eq. (|58D is positive; equality is restored 
when Q is replaced by (5(A). If A > 0, then the left hand side of Eq. (|58D is 
negative. The modification to restore equality is somewhat more complicated 
than just substituting (5(A) for Q and is not given here. 

Hence, from the quantum expectations of the one-body operator (5ij(A), 
its quadratic and cubic scalar couplings, and the value of the nuclear volume 
V , the matrix q and the deformed droplet surface may be constructed. 

5 Nonlinear Collective Model 

The next objective is the construction of quantum collective models corre- 
sponding to the nonlinear gl(3,R) vector fields on Euclidean space. The 
natural solution to the quantization problem is to determine the irreducible 
unitary representations of the relevant real Lie algebra of observables. An 
important requirement is that the Hamiltonian must be a function of the 
algebra observables, and, when this is true, the Lie algebra is called a spec- 
trum generating or dynamical symmetry algebra. For example, the usual 
Hilbert space of square-integrable wave functions that describes a spinless 
nonrelativistic quantum particle is an irreducible unitary representation of 
the Heisenberg algebra generated by the particle position and momentum 
coordinate operators. For a relativistic free particle, the dynamical algebra 
is the Poincare algebra spanned by the Lie algebra of the Lorentz group and 
the four momentum. Both are spectrum generating algebras. 

The spectrum generating algebra for nonlinear collective motion surely 
includes the nonlinear gl(3,R) vector fields. These observables by themselves 
would form a dynamical algebra if the energy were entirely kinetic. But the 
gl(3,R) algebra is insufficient because the potential energy depends upon the 
spatial distribution of particles. Suppose the spatial distribution is defined 
by the nonlinear tensor Qij{A) and its associated deformed droplet. Then 
the collective potential energy may be approximated by the sum of attrac- 
tive surface tension and curvature terms plus repulsive Coulomb energies 
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for a uniform mass and charge distribution whose boundary is the droplet 's 
surface. These potential energies are functions of the parameter A and the 
eigenvalues of the matrix Qij{A) or, equivalently, the rotationally invariant 
functions of the nonlinear Q{A) tensor. This suggests constructing the al- 
gebra generated by the nonhnear gl(3,R) vector fields and the (5(A) tensor. 
Indeed these generators close under commutation to form the general col- 
lective motion algebra gcm(3) and, if exponentiated, the general collective 
motion group GCM(3). 

The general collective motion algebra is generated by two subalgebras. 
The first subalgebra is generated by the Qij{A) tensor and is isomorphic to 
R^, the abelian 6 dimensional Lie algebra. The second algebra is generated by 
the angular and vibrational momenta, L^., Tjj(A), 5'(A), and is isomorphic 
to gl(3,R). The algebra gcm(3)= [R^]gl(3,R) is the span of these two Lie 
algebras and is a fifteen dimensional semidirect sum for which R^ is the 
abelian ideal. Note that two gcm(3) algebras for different A are isomorphic. 

A faithful matrix representation of the general collective motion semidi- 
rect sum is given by the 6x6 real matrices, 

gcm{3)c^l{E,X)^( ^ _f^],E = 'E\, (59) 



where the isomorphism is 



iQij{A) t-> + ej,i+3 

iLk/h ^ eijk{eij + e^+i^s+j) 
iTij{A)/h dj + Cji - e^+j^s+i - e^+i^s+j (60) 

--^(>ij[en + 622 + 633 - 644 " 655 - 6^6) 
i S{A)/?i 1-^ eii + 622 + 633 - 644 - 655 - 666- 

Bkm denotes the 6x6 matrix with 1 at the intersection of row k with column 
I and zero elsewhere. The connected Lie group is given by exponentiation of 
the algebra, 

GCM(3)c^|(A,^)^ (^^ A = *A,yeGL+(3,R)|, (61) 

and obeys the semidirect product multiplication rule 

(Ai, ^i)(A2, 92) = (Ai + g,A2 'gi, gm). (62) 
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Thus, GCM(3) is isomorphic to a semidirect product of the abehan normal 
subgroup of 3 X 3 real symmetric matrices under addition with the general 
linear group GL+(3, R). 

The irreducible unitary GCM(3) representations are given by the Mackey 
inducing construction on associated SO (3) bundles: 

Theorem. For each nonnegative integer C, there exists an irreducible uni- 
tary representation of GCM(3) on the Hilbert space 

i/^ = : G'L+(3,R) ^ C'^+^ | 

ii)-^{gR) = ^{g)V^{R),ior g e GL+{3,R), R e S0{3) 

(ii) [ \\^If(g)\\'du{g)<oc\, (63) 

iGL+ {3,R) J 

where V'-'' denotes the 2C + 1 dimensional unitary irreducible representa- 
tion of S0(3), and du{g) is the invariant measure on GL^{3,Il). Condi- 
tion (i) means that the components of the vector 'if{g) satisfy ^{gR)K = 
J2k'=-c^{9)k'T^k' k{^)- The inner product is 

{<i'{g)mg))diy{g), (64) 

where {v,w) = YS<:=-c '"k'^ k for v,w & 0^*^+^ The action of GL+(3, R) on 
H'-^ is given by 

(7r(x)*)(^) = *(x-V),for x,ge G'L+(3,R). (65) 

The Q tensor acts as a multiplication operator, 

(7r(g,,)*)(^) = {g'g)i^^{g)M g e GL+(3,R). (66) 

Every unitary irreducible representation of GCM(3) is equivalent to a rep- 
resentation for some integral C. Two irreducible representations defined by 
two different integers C are inequivalcnt. 

These irreducible unitary representations model linear and nonlinear col- 
lective motion. But which models are found in nature? The answer to this 
question involves both kinematical and dynamical considerations. First the 
Hilbert space of 74-particle states is a reducible representation of the gcm(3) 
algebra. A model to be kinematically allowed must be an irreducible compo- 
nent in the decomposition of the many-particle representation. Second, for a 
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dynamical symmetry, the Hamiltonian must split into collective and intrinsic 
terms plus vanishing or small cross terms. The collective Hamiltonian acts 
within an irreducible representation while intrinsic terms leave invariant the 
perpendicular subspace to the irreducible representation space. 

The answer to the kinematic problem is that every irreducible repre- 
sentation of gcm(3) is a component in the decomposition of the reducible 
A-particle representation. To prove this assertion a change of variables is 
required for the configuration space R'^'^ from Cartesian coordinates to col- 
lective and intrinsic coordinates. Consider the linear case first (A = 0). The 
general linear group acts naturally on the configuration space and, if yl > 3, 
almost every orbit is diffeomorphic to GL+(3, R). Excluding a set of measure 
zero from the configuration space, a transversal may be chosen that intersects 
each orbit exactly once and whose points (xi, . . . , xa ) e R^^ form a 3^ - 9 
dimensional submanifold M for which 



A 



XnaXnb = ^ab- (67) 



n=l 



Then there is a vector space isomorphism 

jC\GL+{3)) ® C\M) ^ £'(R'^) 



*«)X ^ /, (68) 



where 

f(xi,...,XA) = '^{g)x(xi,...,XA) 



3 



Xr. 



J2 diaXna: (69) 



a=l 



for ^ e jC'^{GL^{3)) and x ^ jC'^{M). This is also a Hilbert space isomor- 
phism, i.e., 

< /1I/2 >=< ^1 ® Xl|^2 ® X2 >=< ^l|^'2 >< Xl\X2 > ■ (70) 

The splitting of the inner product is a consequence of the unitarity of the 
representation of GL^{3, R) on R'^'^ and the invariance of the Haar measure 
on £^(GL+(3, R)). It follows easily that 

U{g){^®x) = {71(9)^) ®X 
Qiji^^x) = {AQij)'^)^X- (71) 
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Thus gcm(3) acts only on C'^{GLj^{?), R)) and leaves the intrinsic wave func- 
tion X alone. Finally the reduction of this representation into irreducibles 
is achieved by making a polar decompostion of the general linear group into 
the product of a symmetric matrix times an orthogonal matrix. A basis 
for £^(G'L+(3, R)) is given by the product of integrable functions on the 
space of symmetric matrices times the representation functions for the or- 
thogonal group Vj^y. Grouping together the 2C -|- 1 basis states for fixed 
C and —C < k' < +C provides the isomorphism with the vector valued 
wavefunctions of the induced representations. Since —C < k < +C is not 
restricted, the multiplicity of each irreducible gcm(3) representation space 
H'~^ in C'^{GL^ (3,11)) equals 2C + 1. Because the space of intrinsic wave 
functions is infinite dimensional, the multiplicity of each irreducible gcm(3) 
representation in £^(R^'^) is countably infinite. 

If A is negative, a similar construction determines the GCM(3) decom- 
position. The principal difference from the linear case is that the separation 
of variables into collective GL^{3, R) and intrinsic transversal coordinates is 
made in the chart space instead of the Euclidean manifold. For A-particles, 
the coordinate system is the Cartesian product of A-copies of Re- 
moving a set of measure zero, is diffeomorphic to the product of GL^{3) 
and the transversal submanifold A/". The Hilbert space isomorphism is 



where denotes the Cartesian product of A-copies of M_ = R^ fl S 



Since the nonlinear action of GL^{3, R) on C'^{M^) is unitary and the mea- 
sure on £^(G'L+(3)) is invariant, the isomorphism satisfies Eq. ([7^). The 
Hilbert space C^{M^) of admissible wave functions is a subspace of £^(R'^'^) 
consisting of the wave functions that vanish for points interior to the sphere. 
Since GCM(3) is represented on the tensor product space in the same way as 
the linear case, Eq. (|7TD , the multiplicity of each irreducible representation 
of gcm(3) in the decomposition of its reducible representation on £"^{^1^) is 
again countable infinite. 

If A is positive, then the necessary separation of variables must be made 
on O^, the Cartesian product of A-copies of the chart space But, since 
the chart space excludes the sphere of radius \/A, the GL^{3, R) orbits in 
are not globally diffeomorphic to the Lie group G'L+(3,R). If i denotes 
the natural injection of into R^"^, then the Hilbert space of states is 
C'^{Q), where the manifold Q = i''^{GL^{3) x A^). Although the regular 



£\GL+{3)) ® C\Af) ^ C\0^) ^ C\M^) 



(72) 
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representation of GL+(3) is not defined anymore because of the sphere's 
obstruction, the Lie algebra gl{3, R) representation, dependent upon the 
local group action, is well-defined on C'^{Q). Since the injection i commutes 
with the derived Lie algebra representation, its decomposion into irreducibles 
is identical to the previous cases. 

For a spectrum generating algebra the collective Hamiltonian must be a 
function of the gcm(3) Lie algebra elements. A phenomenological gcm(3) 
Hamiltonian is provided by a linear combination of rotational scalars formed 
from the generators. For a tractable model the scalars may be restricted 
to polynomials of low degree. Up to quadratic, time-reversal symmetric 
terms, the kinetic energy is a sum of rotational, quadrupole vibrational, 
and breathing oscillation terms constructed from the general linear group 
generators 

T(A) ■ T(A) S{Kf 
^ ~ 21^ 2B + 2C ' ^^"^^ 
where X, B and C are adjustable parameters. A more realistic model would 
include higher degree polynomials including terms like the cubic scalar [T'(A) x 

g(2)(A) X T(A)]0. 

The potential energy is a scalar function of the quadrupole and monopole 
components of the Q(A) tensor. The quadrupole scalars are functions of the 
quadratic [g(2)(A) xg(2)(A)]0 oc /^^ and cubic [(5(2)(a) xQ(2)(A) xg(2)(A)]0 ^ 
cos 87 scalars. Hence, the quadrupole part of the potential may be ex- 
pressed as a function \/(/3,7). The monopole component of the poten- 
tial should have a (5°(A) term. For incompressible motion there may be 
a det(5(A) term too. 



6 Conclusion 

The results of this paper may be amplified, developed and extended in sev- 
eral directions. First, in addition to the GL+(3,R) actions of this article, 
other nonlinear group actions on Euclidean space may be discovered and 
applied to the description of geometrical collective motion in many-particle 
systems. The nonlinear GL+(3,R) collective motion theory developed here 
is an example of the constructive procedure for any nonlinear theory of ge- 
ometrical modes. Although the requirement that the set of collective vector 
fields form a Lie algebra that contains the angular momentum algebra so(3) 
may not be relaxed, a finite dimensional algebra is not essential provided the 
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theory is still tractable, e.g., the collective vector fields might be a Virasoro 
or Kac-Moody algebra. 

Second the adjustable parameters of the collective potential energy may 
be determined from more basic considerations than just phenomenological 
fitting to energy levels and transition rates. For example, the potential energy 
may be defined by attractive surface and curvature energies plus Coulomb 
repulsion for a uniform positively charged fluid contained within the surface 
of the deformed droplets of §4. More fundamentally, the collective energy 
may be derived from a microscopic Hamiltonian H after an intrinsic wave 
function x is chosen. One way to give such wave functions is by projection 
from a Hartree-Fock state '4>hf, 



J (j 



GL+(3,R) 



dn{g) i1^hf{9 ^■xi,...,g ^-xa) 



(74) 



Since /i is the invariant measure on the general linear group, x is an in- 
trinsic wave function, i.e., xi9 ' Xi, . . . , g ■ xa) — x(^i, ■ ■ ■ , xa) for every 
g e G'L+(3,R) and x is well-defined on J\f. The matrix elements of the 
collective Hamiltonian He are defined by 



(75) 



where $, ^ G £^(^^+(3, R)). 

Finally the nonlinear gcm(3) algebra may be extended to the symplectic 
Lie algebra sp(3,R). This algebra consists of the quadratic functions of the 
position and momentum operators or, equivalently, the quadratics in the os- 
cillator phonons. For the nonlinear geometries this algebra should be defined 
using the chart space coordinates of 0±. The position and derivative op- 
erators in the chart coordinates {x'i,X2,x'^) are related to the corresponding 
operators on the Euclidean manifold by 



X, 



d_ 
dx'i 



(1 + 4)^/^ 



d_ 

dxi 



Axi d 

^4 Qj. 



(76) 



Defining the nonlinear phonons by a'^ 
of the real symplectic Lie algebra is 



x[+d/dx'j)/\/2, the complexification 



sp{3, R) = span {af^Ha'^)\ [(a^ta^. + {a'^ya',]/2} . (77) 
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Its maximal compact subalgebra is the nonlinear Elliott u{3) Lie algebra 

t.(3) = span{[(a;)ta; + (a;.)y]/2}. (78) 

The irreducible representatations of the nonlinear symplectic sp{3, R) and 
the nonlinear Elliott su{3) algebra for A-particlcs arc given by highest weight 
arguments similar to the usual linear symplectic and Elliott theories, albeit 
modified to the chart space >C^(C^). These adapted symplectic and Elliott 
models are the natural frameworks for microscopic calculations of nonlinear 
collective motion. 
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List of Figures 

1. Figure 1. Coordinate charts for three-dimensional Euchdean space. 

2. Figure 2. For A = — 1 and a = b = l,c = 2, the droplet is a surface of 
revolution that contains a spherical cavity. 

3. Figure 3. When a = c = 1,6 = 2 and A is positive, the droplet is a 
surface of revolution with a neck. 

4. Figure 4. Triaxial droplets for various positive A and a — l,b — 3,c — 

2. Sections through the principal XZ, YZ, and XY planes are drawn 
for each droplet at the nine, twelve, and three o'clock positions, respec- 
tively. 
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